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t—i ! We consider the effects of "photon bubble" shock trains on the vertical struc- 

! ture of radiation pressure-dominated accretion disks. These density inhomo- 

geneities are expected to develop spontaneously in radiation-dominated accre- 
tion disks where magnetic pressure exceeds gas pressure, even in the presence of 
magnetorotational instability (MRI). They increase the rate at which radiation 
escapes from the disk, and may allow disks to exceed the Eddington limit by a 
substantial factor without blowing themselves apart. To refine our earlier anal- 
ysis of photon bubble transport in accretion disks, we generalize the theory of 
photon bubbles to include the effects of finite optical depths and radiation damp- 
ing. Modifications to the diffusion law at low r tend to "fill in" the low-density 
regions of photon bubbles, while radiation damping inhibits the formation of 
photon bubbles at large radii, small accretion rates, and small heights above the 
equatorial plane. Accretion disks dominated by photon bubble transport may 
reach luminosities of 10— > 100 times the Eddington limit (X_e), depending on 
the mass of the central object, while remaining geometrically thin. However, 
photon bubble-dominated disks with a— viscosity are subject to the same ther- 
mal and viscous instabilities that plague standard radiation pressure-dominated 
disks, suggesting that they may be intrinsically unsteady. Photon bubbles can 
lead to a "core-halo" vertical disk structure. In super-Eddington disks the halo 
forms the base of a wind, which carries away substantial energy and mass, but not 
enough to prevent the luminosity from exceeding Le- Photon bubble-dominated 
disks may have smaller color corrections than standard accretion disks of the 
same luminosity. They remain viable contenders for some ultraluminous X-ray 
sources and may play a role in the rapid growth of supermassive black holes at 
high redshift. 
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Subject headings: accretion: accretion disks — hydrodynamics — instabilities — MHD— 
radiative transfer — X-rays: binaries 

1. Introduction 

Photon bubbles can greatly enhance the escape of radiation from magnetized, radia- 
tion pressure-dominated atmospheres. In photon bubble instability (Arons 1992; Gammie 
1998), the atmosphere spontaneously forms a propagating pattern of low-density channels 
separated by regions of high density. The fastest growing inhomogeneities evolve to nonlin- 
ear amplitude, becoming a train of isothermal, gas pressure-dominated shocks (Begelman 
2001, hereafter B01; Turner et al. 2005). Radiation tends to escape through the underdense 
regions between the shocks, avoiding the regions of high density. The mean radiation flux 
exceeds the local Eddington limit (calculated using the local gravity, g(z): Shaviv 1998) by 
a factor that can approach the ratio of magnetic pressure to mean gas pressure, yet overall 
hydrostatic equilibrium is maintained by magnetic tension and the dynamic cycling of gas 
between high- and low-density phases. Because the shocks tend to evolve toward larger sep- 
arations and higher density contrasts, atmospheres susceptible to photon bubble instability 
can reach very high luminosities. 

Begelman (2002, hereafter B02) argued that photon bubble shock trains may dominate 
the escape of radiation from the inner regions of radiation pressure-dominated accretion 
disks. Such porous disks would be thinner than standard disks of the same luminosity. 
Under extreme conditions, it may be possible for the total luminosity of the disk to exceed 
the global Eddington luminosity by a large factor, without the disk being destroyed by mass 
loss. 

In this paper we attempt to flesh out the effects of photon bubbles on the vertical struc- 
ture of accretion disks, and particularly on the structure of the disks' surface layers and on 
mass loss in the super-Eddington limit. Our analysis is based on analytic and numerical 
calculations of photon bubbles assuming a static background with a uniform magnetic field. 
Further simulations are required before we can answer the crucial question of whether pho- 
ton bubbles form and grow in the presence of magnetorotational instability (MRI), which 
generates a magnetic field that fluctuates in time and space. Arguments based on instability 
growth rates (reviewed in § 3) suggest that the answer is yes. Other issues are susceptible to 
analytic treatment now, and we focus on these. The B02 analysis paid too little attention 
to the effects of finite optical depth and radiation damping on the existence and properties 
of photon bubbles in accretion disks. These effects prove to be crucial. The basic mecha- 
nism of photon bubbles relies on the inverse relationship between radiation flux and density, 
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predicted by the diffusion law. When the optical depth across a region of changing density 
becomes small, the coupling between flux and local density declines. This means that the 
denser regions see a higher flux than predicted by standard diffusion, while the neighboring 
tenuous regions see a smaller flux. Radiation damping (commonly known as "Silk damping" : 
Silk 1968) is a separate effect, and results from radiative diffusion driven by the divergence 
of the velocity field. It is most important at high optical depths. 

The effects of finite optical depth become important when the overall optical depth 
through the atmosphere is still very large, because photon bubbles involve length scales 
much smaller than the scale height of the atmosphere. In § 2 we generalize the nonlinear 
theory of photon bubbles to include the effects of finite optical depth. To do this we require 
a generalization of the radiation diffusion law, to take account of the nonlocal relationship 
between flux and density at low r. The closure we describe in Appendix A is based on 
the model by Ruszkowski & Begelman (2003), modified to include radiation damping. We 
find that the crucial parameter is the optical depth for gas at the mean density, measured 
across a distance equal to the gas pressure scale height. This critical optical depth — tq, in 
our notation — is smaller than the total optical depth by a factor of order the ratio of gas 
pressure to radiation pressure. Radiation damping impedes the existence of photon bubbles 
where the radiation diffusion time through the disk is longer than the gas sound crossing 
time. The ratio of these timescales gives rise to a second parameter, the damping parameter 
£ (B01). (Note that B01 refers to £ as the "radiation drag" parameter. "Radiation damping" 
is a more appropriate term than "radiation drag" in this context, because the underlying 
radiation field is not fixed but rather responds to the flow; see additional comments following 
eq. [7].) Photon bubbles can exist only when £ is sufficiently small; the relationship between 
£ and To determines the photon bubble behavior, as discussed in § 2.2. 

We apply the generalized theory of photon bubbles to radiation pressure-dominated 
accretion disks in § 3. Below a certain accretion rate, corresponding to a luminosity of a 
few percent to tens of percent of the Eddington limit (depending on a and the mass of the 
central object), radiation damping is too strong for photon bubbles to dominate the disk. At 
higher accretion rates, photon bubbles begin to form within the radius where £ ~ 1. Since 
the radiative losses associated with photon bubbles are expected to increase strongly with 
decreasing £ at £ < 1, we argue that there should be a feedback mechanism that keeps the 
disk core at £ ~ 1 at small r. However, we also find that this self-regulated photon bubble 
loss mechanism is subject to the same thermal and viscous instabilities that plague standard 
radiation pressure-dominated a— disks. We reaffirm our earlier result that steady state disks 
can exceed the Eddington limit while remaining thin. 

In § 3.2 we show that disks dominated by photon bubble transport should form a vertical 
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"core-halo" structure. The halo has a much lower density and pressure than the core, yet 
possesses enough optical depth to maintain photon bubbles well outside the core. Where 
the luminosity is sub-Eddington, the core smoothly merges into a homogeneous atmosphere 
at the expected height above the equator. For super-Eddington luminosities the halo forms 
the base of a wind (§ 3.3). We argue that the mass loss, while substantial, is unlikely to 
prevent highly super-Eddington accretion. In effect, the halo insulates the disk core against 
excessive mass loss. In § 3.4, we comment on the effects of photon bubble transport on 
the color temperatures of disks. We summarize our results in § 4, where we also discuss 
important caveats and amplify on two potential applications of the model: to ultraluminous 
X-ray sources and the rapid growth of supermassive black holes. 



2. Photon Bubbles at Finite Optical Depth 

The essence of the photon bubble mechanism is the inverse relationship between radia- 
tion flux and density. At large optical depths the diffusion law for flux is F = — (c//t)Vp/p, 
where k is the opacity (assumed constant) and p is the radiation pressure. In order to main- 
tain hydrodynamical equilibrium, pressure gradients must be relatively weak, but density 
contrasts can be very strong. Thus, radiation flux is excluded from dense regions while 
passing readily through regions of low density. In a photon bubble wave train, the densest 
regions receive so little flux that they are supported against gravity by gas pressure gradients, 
whereas the low-density regions see a super-Eddington flux and accelerate rapidly. 

This picture breaks down where large changes of density occur across regions of low 
optical depth, i.e., where p 2 /t/|Vp| < 1. The flux cannot adjust over a length scale smaller 
than the photon mean free path. Where density gradients are large, the dense regions will 
receive more flux than predicted by the usual diffusion law while the tenuous regions will 
receive less. We stress that we are considering localized density gradients within an otherwise 
optically thick atmosphere. Thus, to a good approximation the radiation field is everywhere 
isotropic. It is the slight anisotropy, giving rise to the flux, that becomes decoupled from 
local changes in the density. 

In this section we consider how the structure of a photon bubble shock train changes 
as the optical depth across the bubble decreases. Effects of finite optical depth set in while 
the total optical depth across the bubble — r A = Pqk\, where p is the mean density and A 
is the spacing between shocks — is still large. We shall see below that the critical optical 
depth, To, is the much smaller optical depth across a gas pressure scale height, r = ponc^/g, 
where c 2 g is the gas thermal speed and g is the local gravitational acceleration. For tq > 1, 
photon bubbles are essentially unaffected by optical depth effects. For r < 1 < r A , the 
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density contrasts within photon bubbles are smaller, but the basic mechanism persists. The 
weakening continues until r A < 1, at which point the photon bubble mechanism fails. 

The remainder of § 2 is more technical than later parts of the paper. Readers who are 
less interested in the physics of photon bubbles may wish to note the definitions in § 2.1.1 
and then skip to § 2.2.3. The principal results needed for the accretion disk discussion in § 3 
are found in Table 1 and § 2.2.4. 

2.1. Equations 

Our analysis is based on a simplified set of radiation-hydrodynamical equations, in which 
the magnetic field is stiff and is inclined to the x— axis by an angle 9, where — n/2 < 9 < n/2. 
We assume that all vectors lie in the x — z plane, and that there is a uniform gravitational 
field — gz. Since the motion is constrained to lie along B, the components of the velocity 
are v = v (cos#, sin#). As in B01, we separate the gas and radiation energetics and assume 
an isothermal equation of state for the gas, with a constant sound speed c g = (pg/p) 1 ^ 2 . For 
additional discussion of the approximations and assumptions behind this model, see B01. 

The five equations that must be satisfied describe continuity, momentum conservation 
parallel to B, conservation of radiation energy, and radiative transfer in both the x and z 
directions. We seek periodic "plane- wave" solutions with wavelength A, in which quantities 
associated with the wave depend on position and time through the combination 

s = x cos + z sin + v p t; (1) 

until further notice we denote differentiation with respect to s by a prime. We may arbitrarily 
assume v p > 0, i.e., that wavefronts move from right to left. Thus, the wavevector makes 
an angle + n with respect to the x— axis, where — n/2 < < ir/2. The quantity v p is 
the phase speed of the wave; note that the wave speed used in B01, v = v p / cos0, is the 
speed of the intersection of a wave front with the x— axis. For purposes of comparison, note 
also the following replacements of the parameters ( and b from B01: ( = tan0, b = tan#. 
According to our conventions, cos0 and cos# are always positive, while sin0 and sin^ may 
take on either sign. 

We consider photon bubbles with wavelengths much smaller than the radiation pressure 
and density scale heights, in which buoyancy forces are unimportant. (See Arons [1992] 
and Begelman [2006] for discussions of photon bubbles driven by buoyancy.) To the required 
level of approximation we may treat the velocity and density as pure functions of s. Defining 
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z-independent wave-averaged quantities by the integral 

(A)^A- 1 [ S+X A(s)ds, (2) 

J s 

we define the mean density (p) = po- In order to maintain overall hydrostatic equilibrium 
of the atmosphere, we cannot ignore the secular variation of pressure with height. For the 
radiation pressure we write p(z, s) = p(s) — pogz, with (p) = p . The "isothermal" sound 
speed associated with the radiation-dominated atmosphere is c = (po/Po) 1 ^ 2 an d H = c^/g 
is the radiation pressure scale height of the atmosphere. 

Continuity . The continuity equation, 

g = -pV-v, (3) 

yields, to lowest order, 

v p p' + cos(9 - (j))(pvy = 0, (4) 

which is readily integrated. To determine the constant of integration, we demand that there 
be no net mass flux through the atmosphere, (pv) = 0. We then have 

p [v p + cos(6> - 4>)v] = p v p . (5) 



Momentum parallel to B . Projecting the momentum equation onto the magnetic field direc- 
tion, we obtain 

OK 

pv' [v p + cos(# — <j))v\ = poVpv' = —c 2 g p cos(# — (/))— pg sin 9 H (F x cos 9 + F z sin 9) . (6) 

Energy . The radiation energy equation can be written 

V • F = -ApV • v - 3^. (7) 

Under conditions likely to apply in accretion disks we may ignore the last term on the right- 
hand side. The first term on the right-hand side leads to damping via radiative diffusion 
("Silk damping" [Silk 1968]; described as the "nondiffusive regime" in section 2.1 of Agol & 
Krolik 1998; see also Blaes & Socrates 2003 for an extensive discussion in relation to linear 
instabilities of radiation-dominated atmospheres). 

Assuming that F is a function of s, we have 



V • F + 4pV ■ v = (F x + Apv cos 9)' cos + (F z + 4pv sin 9)' sin = 0. 



(8) 
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In order to maintain wave-averaged hydrostatic equilibrium in the momentum equation, 
eq. (8) must integrate to 



F z sin + F x cos + 4pv cos(# — 0) 
where Fe = gc/ k is the Eddington flux. 



gc 

K 



sin = F E sin 0, 



(9) 



Radiative transfer . Because we wish to consider the effects of low optical depths within 
photon bubbles, the usual diffusion equation is not applicable. Instead, we will use the closure 
developed by Ruszkowski & Begelman (2003; hereafter, RB), as generalized in Appendix A to 
handle motion. Tests of the stationary (RB) version against Monte Carlo simulations suggest 
that this closure adequately captures the main features of radiative transfer in the presence 
of tilted, plane-parallel density variations deep within an atmosphere. Using equations (A7) 
and (A5), we obtain 



,_, c dp 1 
F r . = -J- + 



pn dx 5pn 



d ( 1 dF x 



+ 



d 



dx \pK dx J dz \pn dx pn dz 



1 dF z 1 dF x 



(10) 



with an analogous equation for F z . Using eq. (8) we can show that 

— 4 cos(# — 0) cos ( 



F x = -^+ 1 



pn dx bpn 
1 



F' 

X 

pn 



pv \ 
pn) _ 



-4cos(#-0)sin0( ^ 

pn 



pn dz hpn i \pn J 
and by using eq. (9) and the adopted form of p(z, s) we have 



V 



g sin + 



Ap KVr 



(Po ~ P)- 



(11) 
(12) 

(13) 



2.1.1. Dimensionless Equations 

At this point, it makes sense to simplify the notation by expressing the equations in 
terms of dimensionless variables. We define: 

V = —; f = — ?; w = — = — ; m p = — — ; r = p KH g ] 14 

po gcsmcp & H g c g \cos(9 - 0) 



R - C l at — c 4^_cos(g-0)j 

P = — ! M = — : 4 = — . (15) 

eg Kp Hc M 



- 8- 



These scalings are physically motivated. We normalize s to the gas pressure scale height 
H g = c^/g, which is the characteristic scale in the high-density region of the photon bubble. 
The quantity m p is a Mach number expressing the speed at which projected wave fronts 
move along a magnetic field line, normalized to the gas pressure sound speed. r is the 
optical depth across a distance H g , assuming that the gas has its mean density po- is 
the diffusivity parameter from Arons (1992), roughly the ratio of the dynamical time to 
the diffusion time across the atmosphere, c/r H c , where t h is the optical depth across the 
radiation pressure scale height. The damping parameter £ ~ thc 9 /c is similar to the quantity 
referred to as the "drag parameter" and denoted by the same symbol in B01; it physically 
represents the ratio of diffusion time to sound crossing time in the gas, measured across the 
radiation pressure scale height. All dimensionless quantities except / are positive-definite. 

In dimensionless form, the radiative transfer equation (11) becomes 



/ = COS0 1 - 



7/ 



1 + 



sin < 



+ 



5r 2 ?? 



£m p cot . 



(16) 



where a prime now denotes differentiation with respect to w. Eliminating F z in favor of F x 
and using eq. (16) to eliminate / in favor of its derivatives, the momentum equation (6) 
becomes 



V 



1\ r , ,„ , N . , tan(0 — 
1 I [cos 0tan(0 — (f>)— t,m p \-\ ^ — 



£m„ cot i 



V 



;i7) 

Equations (16) and (17) describe photon bubbles with finite optical depth and radiation 
damping. 



2.2. Analysis 

Equations (16) and (17) can be combined into a third-order, nonlinear ordinary differen- 
tial equation for r](w). Physical solutions must pass through a critical point at rf = mi and 
a shock transition with the jump condition = m^, where f] + and i]^ are the post-shock 

and pre-shock densities, respectively. The quantities / and f'/rj + £m p cot (p(rj'/r] 3 ) must be 
continuous across the shock front. 

We do not attempt a full numerical solution of equations (16) and (17). Instead, we 
use approximate solutions to estimate the quantitative changes in photon bubble structure 
as t decreases. We restrict our attention to photon bubbles with large density contrasts, 
77_ < K 



- 9- 



2.2.1. Infinite Optical Depth 
In the limit r — > oo, eq. (17) reduces to a single first-order equation in 77, 

(jjT -l)v' = (ri- 1) [cos 0tan(£ - 0) - £] , (18) 

recovering the "stiff-wire" limit discussed in B01 and Turner et al. (2005). 
The critical condition and jump condition imply that 

mj = 77+77- = 1. (19) 

Between the shock and the critical point (where p = po), the gas is supported against the 
effective gravity primarily by gas pressure and the density profile is approximately exponen- 
tial, p ~ p + exp{— u>[cos0tan(6 ) — 0) — £]}. Beyond the critical point, the density varies 
inversely with distance, 77 ps {-u;[cos0tan(# — 0) — The characteristic optical depth in 

the low-density region, t(w) = t i]w = T O [cos0tan(# — 0) — is constant. If t ^> 1, the 
optical depth is large everywhere and we are justified in ignoring finite-r effects. 



2.2.2. Onset of Finite-r Effects 



The effects of finite optical depth appear first near the critical point and in the low- 
density region downstream. Before tackling the case where these effects are dominant, let 
us first consider them in the perturbative limit. To simplify the algebra, in this section we 
consider only the case £ = 0. Let 



77 = 770(1 + £771), / = /o(l + e/i), m 2 p = 1 + 2em 1 , 



(20) 



where e = (5tq) 1 C 1 and 770, fo are the solutions of equations (16) and (17) in the limit 
e — > 0, i.e., 



[l + 770) = -cos0tan(6> - 0), f 



■ COS0. 



(21) 



V6 Vo 

Note that 771 and f\ are functions of w while mi is a constant. Expanding equations (16) 
and (17) to O(e), we obtain 



-if -1) (mvi + vM + 



2(mi - 771)770 



4 



and 



770771 cos + 



770771 COS + 



r)ofo(fi + r)i)- 



tan(0 - 0) (22) 



(23) 
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Equation (22) still has a critical point at 770,0- = 1, /o,cr = 0; from eq. (21) we determine 

that 

1 3 / ' /' V 3 

Vo,cr = cos0tan(# -(f)), rj'^ = - cos 2 0tan 2 (0 - 0), (-M = -- cos 3 0tan 2 (# - 0). 

(24) 

At the critical point the left-hand side of eq. (23) must vanish, while eq. (22) implies mi = 
77 l cr . Therefore, 

3 

mi = ?7i iCr = - cos 2 0tan 2 (# — 0). (25) 
8 

The fact that mi, i]i tCr > means that finite-r effects increase both the wave speed and 
the density at the critical point. These trends occur because extra radiation flux spills over 
into the dense region, compared to the usual diffusion law, when the optical depth is finite. 
Consequently, radiative acceleration takes over from gas pressure support at a higher density, 
pushing the flow through the critical point sooner and propelling the wave. 

We next consider the flow well downstream of the critical point. For 770 and fo we use the 
asymptotic solution discussed in § 2.2.1: r] = [w cos0tan(# — 0)] _1 , fo = — w cos0tan(# — 0). 
Assuming 770 <C 1, equations (22) and (23) yield wq^ + rji + \ cos 2 0tan 2 (# — <fi) = 0. Given 
that rji is pinned to ?7 1)Cr at w ~ a few, the solution at w ^> 1 must converge to a constant 
value, 

77i = -icos 2 0tan 2 (# - 0). (26) 

Thus, at large but finite r , the density is depressed below its infinite-r value at w ^> 1. 
However, we will see below that this behavior reverses as tq decreases to the point where the 
finite-r corrections are no longer perturbative. At small To, 77 oc (tqw)^ 1 becomes very much 
larger than its infinite-r value. 

Although the density at w ^> 1 is decreased by finite-r effects, the minimum density 
of the photon bubble, is higher, given a fixed maximum density 77+. The shock jump 
condition t?_ = m 2 /7] + implies that increases with m 2 Thus, finite-r effects tend to wash 
out the density contrasts of photon bubbles, a trend we shall continue to find as To becomes 
small. 



2.2.3. Arbitrary Optical Depth 

Guided by our results above, which show that the 77 oc w^ 1 scaling is preserved to 
lowest order when finite-r effects are included, we seek a solution of the form 77 w A/w well 
downstream of the critical point, i.e., when 77 m p . In this limit, both (f'/r))' and (rj'/rj 3 )' 
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are constant. To model solutions that pass through the critical point, we extrapolate the 
assumed A/w dependence to the critical pont at r\ = m p , and demand that the right-hand 
side of eq. (17) vanish at this point. This implies that 



1 - 



tan(0 - 0) 



5tq [cos0tan(6> — 0) — £m p 



+ £m p cot i 



7[ 



(27) 



Extracting the factor 1 — m p /rj from eq. (17) and returning to the limit i] <C m p , we are left 
with 



— Tfirp—— = cos0tan(# — (f))— t;m v . 
r] 2 

Setting r] = A/w, eq. (28) allows us to express A in terms of m p , 



(28) 



cos tan(0 — (f))— £m p 



(29) 



while equation (16) gives 

f' = 



sin 9 



£m p 



m; 



sin0sin(6> - 0) A A 2 tan(6»-0)' 



(30) 



Since A must be positive, solutions are possible only for cos0tan(# — 0) > C,m p > 0. 
Substituting from eq. (30) into eq. (27), we obtain 



m. 



1 + 



m r . 



5t$A 2 



{ 1 + ^m~ i 1 + 2cot ^ tan (^ -0)]} 



(31) 



Equation (29) can then be used to turn this into a quadratic equation for m p . 

If £ min[cos 0tan(# — 0), V5t ], we may ignore the effects of damping (i.e., set £ = 0). 
Equation (31) then has the solution 



m r 



1+1 + 



4cos 2 0tan 2 (6' - 0) 
5r| 



1/2" 



(32) 



In the perturbative limit r ^> 1, we obtain a value for mi greater than eq. (25) by a factor 
8/3, but the angular dependence is the same. The discrepancy in the coefficient reflects the 
crudeness of our method for approximating solutions that pass through the critical point. 

To simplify the analysis further, we specialize to the case 9 = 0, which is appropriate 
for conditions inside accretion disks where the predominant magnetic field is expected to be 
roughly parallel to the equatorial plane (Balbus & Hawley 1998). We then require sin0 < 0, 
so that tan(0 — 0) > 0. (This restriction implies that photon bubble wavefronts propagate 
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upwards, as seen in simulations [Turner et al. 2005].) Defining ji = |sin0| and using 
e = (5tq )~ 1 from § 2.2.2 without assuming its magnitude, we obtain the quadratic 

m\{\ - ei 2 ) + m p [2e^ - (1 + £ £ 2 )] - £//(// - = 0. (33) 

Since we identify the critical point with the effective sound speed (which can be boosted 
above the gas sound speed by radiation pressure), we seek solutions with 1 < m p < 
where the upper limit guarantees that A > 0. This means that there are no photon bubbles, 
or very weak ones, for £ > fi — the radiation damping is too strong. We note that the 
same restriction applies to the development of the linear photon bubble instability (Gammie 
1998; Blaes & Socrates 2001, 2003). Assuming that /i ~ 0(1), we may solve eq. (33) 
approximately for the three characteristic orderings of r , £, and 1 that admit nonlinear 
photon bubble solutions according to our simple ansatz. The properties of these solutions 
are shown in Table 1. 

Case (a), labeled "ideal" , was treated in B01 and Turner et al. (2005). In this solution the 
standard diffusion approximation is valid everywhere. The other cases involve modifications 
to the diffusion law, due to low optical depths. Case (b) illustrates the effects of low r in the 
limit where radiation damping is unimportant. In this case m p ^> 1, continuing the trend 
toward larger wave speeds that we saw in § 2.2.2. A higher wave speed, in turn, implies that 
the density contrast is smaller, from the shock jump condition = mi with (rj) = 1. 

Also, note that A « (VStq)^ 1 ^> 1, reversing the trend toward lower density that we found 
in the perturbative limit. The effects of low optical depth cause the low-density region of 
the photon bubble to "fill in", at least in the undamped limit. Note that the characteristic 
optical depth across a density scale length is At ~ 5 -1 / 2 . Thus, the optical depth across 
the low-density part of the wave pattern is regulated to be ~ 1. Case (c) represents another 
situation in which the wave speed can be much larger than the gas sound speed. In this case, 
radiation damping approximately balances the accelerating forces, hence the label "critically 
damped" . A is large, and the bubble is once again filled in — here the characteristic optical 
depth is at least ~ 1 and may be much larger. In the limit m p approaches 1 but the 



Table 1. Parameter Space for Photon Bubbles 



Case 


Ordering 


m p 


A 


Description 


a 


(«1«T 


1 


l//i 


Ideal, undamped 


b 


£ < TO < 1 




1/V570 


Diffusion-dominated, undamped 


c 


TO < £ < 1 




max [l/v^ro, /i/£ 2 ] 


Diffusive, critically damped 
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density remains large provided that r <C £. 

We can use the expressions for m p and A to estimate the wavelength of the photon 
bubble as a function of the density on either side of the shock front. Assuming that we may 
extrapolate the density profile rj « A/w all the way to the shock, we have A w A/r}-. The 
shock jump condition then gives 

A » -^rj + . (34) 

It will prove convenient to express A in terms of the maximum density Simulations sug- 
gest that nonlinear photon bubbles tend to evolve toward the longest permitted wavelength 
(Turner et al. 2005), corresponding to the largest permitted value of r] + . The latter is either 
set by the magnetic field strength, which must be large enough to resist buckling under the 
postshock gas pressure, or the wavelength, which must be smaller than the scale height. 



2.2.4. Mean Flux 



For the purpose of studying the release of radiation by accretion disks, we are interested 
in the vertical flux as measured in the laboratory frame, 



(i^iab) = (F z + Apv sin 9). 



(35) 



Since the local flux increases linearly with w in the low density part of the bubble, the mean 
flux is given approximately by the local flux at w = A/2. Since we are assuming 9 = 0, the 
second term on the right-hand side of eq. (35) does not contribute. From equations (17) and 
(12), we obtain the Eddington enhancement factor 
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3. Application to Accretion Disks 

Can nonlinear photon bubbles form in accretion disks, and if so, do they affect the 
disk structure and luminosity? The answer depends largely on the relationship between the 
dissipation rate and the density as a function of height above the disk equator. Simulations 
of magnetohydrodynamic disks, both with (Turner, Stone, & Sano 2002; Turner et al. 2003; 
Turner 2004) and without (Miller & Stone 2000; Hirose et al. 2004) radiation transport, sug- 
gest that most of the dissipation may occur in a magnetically dominated corona, extending 
to several density scale heights. If most of the radiation is produced so high in the corona 
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that the optical depth is low (e.g., Socrates & Davis 2005, and references therein), then the 
role of photon bubble transport is moot. If, however, the bulk of the energy is produced 
sufficiently close to the equator that r is large, then photon bubbles can play a key role 
in transporting radiation to the disk surface. The least favorable case for photon bubble 
transport is the one-zone slab model, in which most of the dissipation occurs in the dense 
core of the disk. In this paper we restrict our attention to the latter model. B02 considered 
this case but did not take into account the effects of radiation damping, which prove to be 
important at large radii and close to the equator. In § 3.1, we will re-analyze the effects of 
photon bubbles on the structure of disk cores, in the light of our analysis in § 2. 

Crudely, photon bubble shock trains can be treated as localized phenomena within 
radiation pressure-supported atmospheres. Since they propagate at roughly the gas sound 
speed, they would require many orbital times to cross the disk scale height. MRI rearranges 
the magnetic field on a timescale of order the orbital time, suggesting that photon bubble 
wave trains are continuously being destroyed and reformed before they have time to traverse 
great distances. On the other hand, photon bubbles typically grow to nonlinear amplitudes in 
less than an orbital time (see below). Simulations show that short wavelength modes quickly 
evolve toward longer wavelengths and larger density contrasts, until growth is quenched by 
either buckling of the magnetic field or lack of space (i.e., when the wavelength exceeds the 
scale height) (Turner et al. 2005). Where shock fronts extend far enough in height to sample 
a range of disk conditions, they curve in order to track the local wavevector of maximum 
instability. These properties allow us to consider the effects of photon bubbles on disks as 
a function of height z above the equatorial plane. In § 3.2 we will consider the gradual 
disappearance of photon bubbles with height, as the atmosphere becomes too tenuous to 
support large local variations in radiation flux. Even if the core supports large Eddington 
factors I S> 1 due to the presence of photon bubbles, above some height the disk must behave 
like a homogeneous disk — material is levitated until the local flux equals the Eddington 
flux in the effective gravitational field. The transition from a porous to a smooth atmosphere 
is especially interesting if the core flux exceeds the Eddington limit: in this case, material 
from the upper layers of the disk cannot remain in equilibrium and is blown off in a wind. 
In § 3.3 we will sketch the possible nature of this wind, and the extent to which it is likely 
to disrupt the accretion flow. 

Because the magnetorotational instability (MRI) leads to turbulence, with magnetic 
field fluctuations on timescales of order the orbital time, there is uncertainty as to whether 
photon bubbles can actually grow and saturate in accretion disks. Several authors have 
argued that photon bubbles are likely to develop even in the presence of the time-varying 
field, because the linear growth rates of photon bubble instability are typically faster than 
MRI growth rates (Gammie 1998; Blaes & Socrates 2001, 2003; B01; B02; Turner et al. 2005). 
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To quantify this statement, we consider the maximum photon bubble growth rate using Blaes 
& Socrates (2003, eq. [93]), in the limit of tight thermal coupling, a stiff magnetic field, and 
electron scattering opacity. In terms of our variables, the imaginary part of the frequency is 

u, = ^ |C °y~ 0)l [cos 0| tan(0 - 0) | - £] . (37) 

ICg 

This rate is an asymptotic limit for wavelengths A < H g ; at longer wavelengths the growth 
rate declines oc A -1 / 2 . For the interior of an accretion disk with negligible damping, the 
maximum growth rate is roughly (3~ l l 2 VL/2 (for a vertical field), which is much larger than 
the maximum MRI growth rate 3fi/4 in the presence of a vertical field (Balbus & Hawley 
1998). For a horizontal field the maximum photon bubble growth rate is a factor of 2 smaller. 
The maximum photon bubble growth rate only becomes comparable to the MRI growth rate 
for wavelengths approaching the disk scale height. 

It is not clear whether long-wavelength (A S> H g ) photon bubble shock trains — those 
associated with large density contrasts — grow directly from linear instabilities, or from the 
merging of shorter wavelength shock trains. Simulations (Turner et al. 2005) seem to indicate 
that the latter process might dominate, in which case long wavelength shock trains could 
develop even more quickly than predicted by the extrapolation of linear stability theory. 

We note that the numerical growth rates presented in Turner et al. (2005: especially 
Figures 8, 11, and 13) are normalized to Q, and give the impression that photon bubble 
growth rates do not exceed a few times fl For purposes of comparison with MRI, this is 
misleading, because the results presented correspond to a case with rather strong damping, 
and therefore strong suppression of the growth rate. Across the disk surface layer considered 
by Turner et al. (2005), £ (= 4| cos(# — <f))\c g E/3F in the notation of Blaes & Socrates 
2003 and Fig. 3 of Turner et al. 2005) decreases with height from 0.53 to 0.43 for the modes 
plotted in Fig. 8, which have = 0, = 63°. This is close to the fastest growing mode under 
these conditions, and the growth rate is ~ (1.5 — 2)0. If damping were negligible, the fastest 
growing mode for 8 = would have = 45°, and the growth rate would be ~ 5.2f2. Thus, 
damping cuts the growth rates in the Turner et al. (2005) simulation by a factor ~ 2 — 3. 

In addition to time-dependence, MRI and other instabilities create fluctuations in mag- 
netic field direction and strength. Photon bubbles are insensitive to magnetic field strength 
provided that the magnetic pressure exceeds the maximum gas pressure in the wave train. 
But fluctuations in direction will distort and refract photon bubble wavefronts (depending 
on the fluctuation scale compared to the wavelength), and may lead to large fluctuations in 
photon bubble growth rate and disk porosity. Where damping is present, a sudden change in 
field direction could damp out (or enhance) a wave train. These effects are beyond the cur- 
rent state of analytical models. Effects due to fluctuations in field direction may be largely 
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mitigated by the fact that MRI leads to a predominantly azimuthal field (e.g., Turner 2004), 
with relatively modest fluctuations in other directions. Moreover, a strong azimuthal field 
can slow (Blaes & Socrates 2001) or even quench (Pessah & Psaltis 2005) the growth of MRI 
modes. 

In the absence of further guidance from simulations, we will assume that localized photon 
bubble wave trains grow wherever possible given the local conditions, and saturate at close to 
the maximum permitted amplitude. All other factors being equal, the amplitude is limited 
by the magnetic field, which must be strong enough to resist buckling. The maximum 
dimensionless density is i] + ~ p m /pg, where p m is the magnetic pressure and p g = poc 2 is 
the mean gas pressure (B01; Turner et al. 2005). In our disk models we will find that 
other factors may limit the amplitude of the photon bubbles, such as constraints on their 
wavelengths (which must be shorter than the scale height) or the transition of the disk to a 
parameter regime in which photon bubble formation is inhibited. Reassuringly, our models 
for disk cores will indicate such a high efficiency of radiative transport (where damping is 
unimportant) that saturation at the maximum permitted levels is unnecessary. 

In the models presented below we will take 9 « 0, corresponding to a field parallel 
to the disk equator. Simulations suggest that (B 2 ) ~ 3(5^,5 r ) (e.g., Hawley 2000), so we 
will adopt p m ~ 3apo, where a is the usual viscosity parameter and po is the radiation 
pressure. Simulations in an initially static, nearly horizontal magnetic field (Turner et al. 
2005) suggest that the dominant shock trains have the same orientation as the fastest- 
growing linear mode (Blaes & Socrates 2003), corresponding to ~ —70° in our notation 
when radiation damping is unimportant. We will adopt this value for 0, except where 
specified; note that the dominant modes are likely to be different when damping and/or 
finite-r effects are important. 

3.1. Photon Bubbles in Accretion Disk Cores 

Following B02, we adopt a one-zone model in which most of the disk luminosity is 
produced within the high-density core of the disk. We note, however, that the scaling of 
the damping parameter with height (assuming LTE) is £ oc p 9 ^ 8 /z, implying that £ is very 
large near the equator and that photon bubbles are likely to be inhibited in this region. 
The dependence of £ on the orientation of the photon bubble wave train (£ oc cos(0 — 0) ~ 
(1 — yU 2 ) 1 / 2 ; eq. [15]) could partially compensate for the increasing effects of damping in the 
interior of the core. Here we will assume that most of the energy is generated in the upper 
half of the core, and adopt \i = sin 70°. 
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B02 equations (5) through (12) are unchanged. B02 equations (2) and (3), which repre- 
sent an early attempt to model the undamped cases (a) and (b), must be modified: the term 
p 2 CgK 2 /g 2 in (2) and p 2 K 2 / g 2 in (3), which equal Tq , must be replaced by a single power of t . 
The correct formula has one power of To, not two, because the wave Mach number m p oc r^ 1 
in the low-r limit, according to our analysis in § 2. B02 incorrectly assumed that the wave 
speed was unchanged by finite-r effects. In any case, we shall show below that cases (a) and 
(b) occupy only a small region of the accretion disk parameter space — therefore, we shall 
have to replace the entire B02 analysis by a revised treatment based on § 2 of this paper. 

As in B02, we assume a steady accretion disk and normalize the accretion rate to the 
Eddington rate, m = M/Me, where Me = 4hGM/k,c and M = mM Q is the mass of the 
central object. We normalize r to the gravitational radius, x = rc 2 /GM, set a = a/0. 03, 
and define V = 1 — (x- m /x) l l 2 for the standard assumption of zero torque at the innermost 
stable orbit, x- m (Shakura & Sunyaev 1973). We eliminate the dimensionless disk thickness, 
5 = H/r, in favor of the Eddington factor £, by using B02 eq. (7), £ = 3mD/8x. A standard 
accretion disk, without photon bubbles, has £ — 1. The effect of photon bubbles is to increase 
£ so that the disk can lose more radiation without thickening. 

As we showed in § 2, the key parameters determining the properties of photon bubbles 
are the characteristic optical depth across a gas pressure scale height, tq = ponrSfi, and the 
damping parameter, £ = 4/3 1//2 cos(# — (ft) /Mo ~ A(3 1 ^ 2 £ cos(# — (ft) /3a. Using the formulae 
from B02 we have 

r w 4.3 x lO-'a-^imV^m' 1 ^ 17 ^ 33 / 8 , (38) 
^3.6x 10- 2 a- 9 /VP) _1 ™~ 1/8 * 17/8 z 21/16 . (39) 

Figure 1 shows the parameter space of disk core models with £ — a — 1, X[ n — 6, and two 
values of the central mass. The upper excluded region corresponds to 5 > 1, i.e., disks that 
are so thick that they must lose mass. Since this figure assumes £ = 1, it corresponds to 
standard disk models without photon bubbles — we can use it to determine when photon 
bubbles may form. The lowest region corresponds to £ > 1, implying that photon bubbles 
are inhibited by radiation damping. The regimes assumed by B02, the undamped cases 
(a) and (b), occur only in the narrow strips between the curves. Most of the parameter 
space that supports photon bubbles lies between the upper curve and the excluded region, 
and corresponds to case (c) - - critically damped bubbles. Photon bubbles can exist for 
parameters such that £ < 1 and 5 < 1. For fixed values of m and a, the size of this region 
in the m — x plane is a rapidly decreasing function of £. We can see this in Fig. 2, which 
shows how quickly the parameter space available for photon bubbles becomes restricted with 
increasing £. At a given x, all boundary curves in the parameter plane move toward higher m, 
but not at the same rate. The 5=1 curve and the r = £ curve rise approximately linearly 
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Fig. I. — Parameter space of thin (5 < 1), standard (£ = 1) disk cores with a = 1 and 
M = 10M Q (left panel); M = 1O 8 M (right panel). The upper shaded region in each plot 
violates the thinness condition. The lowest (thickest) curve in each plot corresponds to 
£ = 1, the next curve corresponds to ro = 1, and the highest (thinnest) curve corresponds 
to £ = ro. Letters on the right-hand panel refer to the photon bubble parameter space 
discussed in § 2 and summarized in Table 1. Photon bubbles cannot exist in the region 
below the lowest curve. Between the upper curve and the excluded region, photon bubbles 
exist with radiation damping approximately balancing driving forces. The analysis by B02 
is valid only in the narrow strips a and b between curves. 

with £, while the £ = 1 and r = 1 curves increase roughly as rh oc £ 2 and oc £ 3 ^ 2 , respectively. 
The maximum possible £ is therefore reached when the disk conditions approach £ ~ 1, i.e., 
close to the value of £ where photon bubble formation is cut off. This implies 

£ < 4>a* ~ 5a 9 / 17 {mVf^m^x' 21 ^. (40) 

Photon bubbles can affect the disk only where £ max > 1, which implies that photon bubbles 
are important only for 

x<x ph ^13a 18 ^m w / 21 m 2 / 2 \ (41) 

where we have assumed T> w 1. At the same time, the mass flux that can be carried by the 
disk is constrained by the thinness condition 5 < 1, which implies 

£ > £min = l.hmVx- 1 . (42) 

Self-consistency demands £ min < £ max , and therefore that 

rh < Vam l ' 9 V-V l \ (43) 

We can place these constraints into a consistent picture for the role of photon bubbles 
in the central regions of an accretion disk. For rh < S^a^^m' 1 / 8 , photon bubbles never 
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Fig. 2. — Parameter space of thin (5 < 1), disk cores with a — 1, £ — 10 and M = 10M© 
(left panel); M = 1O 8 M (right panel). The upper shaded region in each plot violates the 
thinness condition. The dashed line shows the thinness constraint for £ — 1. The lowest 
(thinnest) solid curve in each plot corresponds to £ = tq > 1; the intermediate thickness 
curve corresponds to tq — 1. In the right panel, photon bubbles are possible only between 
the heavy curve (£ = 1) and the shaded region (case c). In the left panel the heavy curve 
has moved into the shaded region, implying that disks with M = 1OM , £ = 10 violate the 
thinness constraint. Note that all boundaries move upward with increasing £, thus restricting 
the parameter space available for photon bubbles. 

form in the disk — £ > 1 at all radii. For an accretion efficiency of e ~ 0.1, this limit 
corresponds to disks less luminous than ~ 0.3a~ 9 ^ 8 m~ 1 ^ 8 LE- Note that photon bubbles set 
in at lower luminosities (relative to the Eddington limit) for higher mass black holes. If 
m > 3.2<5: -9 / 8 m _1 / 8 , then photon bubbles start to become important at x < rc p b, the radius 
at which the disk first "hits" the £ = 1 curve. Since nonlinear photon bubbles tend to 
evolve toward larger density contrasts, and therefore higher £, a photon bubble-dominated 
disk core should have a natural tendency to track along the boundary of marginal stability, 
corresponding to the curve £ = 1. In other words, the tendency for £ to be as large as 
possible supplies feedback that drives the disk toward £ ~ 1. Since the location of the £ = 1 
curve in the m — x plane depends on £, the value of £ increases with decreasing x, lifting the 
£ = 1 curve and allowing m to remain constant. 

The Eddington enhancement factor ~ 1 at x p ^ and increases as £ ~ (x/x ph )' 21 ^ 4 toward 
smaller radii, provided that the disk is able to sustain the required density contrasts. The 
geometric disk thickness 5 oc (rrf)" 1 oc x~ 13 ^ u at fixed m (and V xs 1); therefore the disk 
thickens with decreasing radius. 

If rh < 170am 1 / 9 the photon bubble-dominated disk does not violate the thinness cri- 
terion at any radius. Therefore, it may be possible for photon bubbles to liberate radiation 
at a rate > 20 — 130(e/0.1)aL E , for black hole masses ranging from 10 — 1O 8 M . Of, course, 
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such supercritical luminosities would drive off mass from the outer layers of the disk where 
the photon bubble mechanism fails — we will address the consequences of this mass loss in 
§3.3. 

For higher accretion rates, £ increases with decreasing x as ~ (x/x ph )~ 21 ^ 34 , until S ~ 1. 
At smaller radii the disk loses mass at a rate sufficient to keep S ~ £ ~ 1, i.e., m oc a; 13 / 18 
where T> 1. Note that the gravitational energy release still increases with decreasing 
radius, despite the mass loss. The decline of rh continues until x 13 ^ 18 !)^ 1 is minimized, at 
x ^ 17. Within that radius the accretion rate levels off at the maximum value calculated 
above, m ~ 1705m 1 / 9 . 

The foregoing discussion assumes that photon bubbles grow to whatever density contrast 
is needed in order to keep the disk near the threshold £ ~ 1 and £ ~ l masL - Simulations 
(Turner et al. 2005) suggest that density contrasts in nonlinear photon bubbles tend toward 
the maximum values permitted by either the magnetic field strength or the wavelength 
(which cannot exceed the scale height). For all cases with £ > 3 — 6 (where the threshold 
depends on a and m), we are in photon bubble case (c), with r C 1 and £ ~ fi. Equation 
(36) then gives 

u 2 

£ ~ 1 + — . (44) 
2^_ 

For reasons mentioned earlier we adopt <fi = —70°, so fi = 0.94 and £ ~ 0.4/ry_ in the photon 
bubble-dominated regions of the disk. Since m p « 1 and = m 2 from the shock jump 

condition, the maximum density in the photon bubbles must be i] + ~ 2.3£ times the mean 
density. 

As discussed in the introduction of § 3, the maximum density is limited by the stiffness 
of the magnetic field to i] + < 3a/ [3, where f3 is the ratio of mean gas pressure to radiation 
pressure. Is this limit on i] + consistent with the value required to maintain £ = 1 and 
£ = £ max ? For x in = 6, half the accretion disk luminosity (for fixed fn) is produced outside 
x\/2 = 23.5, and we will adopt this as the fiducial radius for the purpose of this argument. 
At X1/2 the maximum accretion rate consistent with the thinness criterion is m ~ 180am 1 / 9 , 
corresponding to £ ~ 5.6am 1 / 9 and r] + ~ 13am 1 / 9 . Substituting these values of m and £ into 
f3(x 1 / 2 ) ~ 0.1a~ 1 / 4 m~ 1 / 4 m~ 2 £ 9 / 4 , we obtain the constraint r] + < 600am 2 / 9 , which is easily 
satisfied. Thus, the strength of the magnetic field is more than adequate to support the 
required level of density inhomogeneity. 

The other relevant constraint is that the characteristic separation of the dense regions 
(the photon bubble wavelength) be shorter than the disk scale height, H. From eq. (34), we 
have \/H ~ (3i]+/ V^tq. For the maximum accretion rate at X1/2 we have tq ~ 0.044m -1 / 9 , 
implying X/H ~ 0.02a. The dominant scale length of photon bubbles is therefore much 



- 21 - 



smaller than the pressure scale height. 

3.1.1. Thermal and Viscous Stability 

Domination of radiative transport by photon bubbles does not appear to cure the ther- 
mal and viscous instabilities that plague standard a-disks when radiation pressure exceeds 
gas pressure (Lightman & Eardley 1974; Shakura & Sunyaev 1976). Consider a radiation- 
dominated disk with surface density £ and pressure scale height H . In hydrostatic equilib- 
rium the radiation pressure scales as p oc Eif while the density obeys p oc £/if and the gas 
pressure (in LTE) satisfies p g oc £ 5 / 4 /# 3 / 4 . Therefore, f3 oc E 1 / 4 ///?/ 4 and f oc E 9 / 8 ^ 1 ^ 
The rate of dissipation per unit area, according to the a-model, scales as Q + oc if 2 £, while 
the rate of energy loss scales as Q~ oc £H. 

Consider thermal stability first. Since thermal timescales are much shorter than viscous 
timescales, we can consider E to be constant. In the absence of photon bubbles, £ = 1 and 
Q~ increases with H more slowly than Q + , leading to thermal instability (see, e.g., Frank, 
King & Raine 2002). When photon bubbles dominate the energy loss, however, the disk 
must lie close to marginal stability, £ ~ /i; otherwise the natural tendency for £ to be as 
large as possible would release too much radiation. At constant £, £ ~ r H c g /c oc H l l s in 
LTE. Therefore, a slight increase in H presumably leads to a significant decrease in £ as 
£ increases, thus making the disk less susceptible to photon bubble instability. Since Q + 
increases oc H 2 while Q- probably decreases (or at least increases more slowly), the disk is 
thermally unstable. Note, however, that the dependence of £ on H is very weak, and depends 
entirely on the slight increase of gas temperature with radiation pressure under LTE. The 
conclusion would be reversed if £ decreased slightly with increasing H, either through a 
decrease in gas temperature (unlikely), a small decrease in £ (since £ is much more sensitive 
to S than to H), or a decrease in the characteristic value of /i. 

Now consider viscous instability. According to the a-model, the viscous couple satisfies 
G oc pH oc H 2 H. On viscous timescales, thermal equilibrium is maintained, implying 
Q + = Q~ , where Q + oc H 2 lu. When photon bubbles are absent, Q~ oc H and thermal 
equilibrium gives H oc E -1 , implying G oc E _1 . Since the viscous couple is a decreasing 
function of surface density, the coefficient of viscous diffusion is effectively negative and the 
disk is unstable to clumping into rings (Lightman & Eardley 1974). When photon bubbles 
dominate, £ ~ constant. The viscous couple decreases rapidly with E, and the disk is 
unstable. 

Photon bubble transport could quench both thermal and viscous instability if the density 
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contrasts were close to their maximum values, or scaled with these values. In this case, 
£ oc p/pg oc p 3//4 /p oc iJ 7 / 4 ^" 1 / 4 in LTE (the LTE assumption is not crucial to the qualitative 
result). We therefore have Q~ oc iJ n / 4 X -1 / 4 , cooling increases more rapidly than heating 
at fixed S, and thermal instability is quenched. By the same token, in thermal equilibrium 
the viscous couple increases strongly with S, G oc X 13 / 3 , and the disk is extremely stable on 
viscous timescales. 



3.2. Disk Halo 

We characterize the disk halo as the plane-parallel region, sandwiching the core, through 
which a constant mean radiation flux passes. We assume that photon bubbles also develop 
within the halo. In LTE, we have the scalings £ oc p 9 ^ 8 /z, tq oc pp 1//4 /z, and (3 oc p/p 3 ^ 
(qualitatively similar to the isothermal scalings £ oc p/z, tq oc p/z, and (3 oc p/p). 

Let us denote core quantities by a subscript "c" and use a hat to label quantities 
normalized to core conditions. Since the flux is constant while the Eddington flux increases 
linearly with z, the condition of hydrostatic equilibrium requires that £ = £/£ c = {z/H c )~ l = 
z-\ 

Although £ c ~ 1, £ decreases rapidly with height above the core. The halo departs from 
marginal stability and therefore photon bubbles should be able to approach the maximum 
permitted density contrast. As the atmosphere becomes increasingly porous, both the density 
and pressure drop steeply. This transition zone ends where the maximum allowed value of 
£ approaches £ c . Since the density drops at least as steeply as the pressure, the ordering 
t <C £ is preserved and we remain in case (c) of the photon bubble parameter space with 
A ~ l/v^To. The Eddington enhancement factor is still given by eq. (44). As discussed 
earlier, the magnetic field strength places an absolute upper bound on the postshock density 
of 7] + < 3a/ /3, and therefore a lower bound on However, this upper bound is superseded 
by the constraint that the photon bubble wavelength, A ~ A/rj- (in units of the gas pressure 
scale height), not exceed the atmospheric scale length (or a more stringent upper limit 
associated with coherence length of the magnetic field). Setting the maximum wavelength 
to 0.3xz, where x < 1 is a constant parameter intended to mimic the effects of an incoherent 
magnetic field, we obtain £ ~ 0.3xP KZ , implying 

^~ 10 X _1 ^" 2 «1 (45) 
in the region above the transition zone. 
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In scaled variables the equation of vertical hydrostatic equilibrium becomes 

dp _ | 

— = — zp oc z (46) 

dz 

at z <C r, implying that p and p drop by roughly the same factor in the transition zone, 
but that p decreases with height logarithmically at higher z. This behavior continues to the 
"top" of the halo at 

z max ~ min [4, 5~ l ] , (47) 

where S c = H c /r is the geometric thickness of the disk core. If £ c < then the flux 
passing through the halo at x is sub-Eddington with respect to the mass of the black hole. 
Zmax is then the height at which the optical depth across the halo, p(z max )Kz mSuX , declines 
to ~ 0(1), and photon bubbles can no longer be supported. z max is also the thickness of 
an "equivalent" standard accretion disk with the same local flux but no photon bubbles. 
The difference between the disk-halo structure dominated by photon bubble transport and a 
standard accretion disk is that most of the column density in the photon bubble-dominated 
disk is found in the much thinner core — only a fraction £~ l of the disk column is found in 
the upper regions of the halo. Nevertheless, the scattering photosphere is located near the 
top of the halo, and therefore the disk will appear to be as thick geometrically as a standard 
disk of the same luminosity. 

If 5" 1 < £ c , then the halo flux is super-Eddington with respect to the black hole mass, 
and z max > r. In this case the hydrostatic disk halo is replaced by a radiation pressure-driven 
wind. We discuss the properties of this wind below. 



3.3. Disk Wind 

If the disk core thickness H c exceeds £~ 1 r (i.e., for £ C 5 C = 3mV/2x), the disk is super- 
Eddington and must lose mass. This condition is equivalent to the halo being geometrically 
thick, z max > r. In this section we estimate the properties of the wind and show that the 
mass loss need not be catastrophic. 

Treating the outflow is more complicated than treating the static halo because part of 
the energy flux is advected, F a d v = 4pi> (see § 2.2.4). Naively, this advective flux seems 
likely to dominate the energy budget in the wind, leading to such a large mass flux that it 
would decimate the accretion disk. To see this, suppose that the halo carries a conserved, 
vertical wind flux fn w = pv w at z < r. At z ~ r, the gravitational acceleration switches from 
increasing with height oc z to decreasing oc r~ 2 . Let us assume that the flow reaches escape 
speed near this point, v w ~ vk — (GM/r) 1 / 2 , where vk is the Keplerian speed at r. We 
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saw in § 3.2 that the static halo (above the transition layer) has roughly uniform radiation 
pressure, p ~ l^x^Pc^cX' 1 ^ 2 (to within a logarithmic factor). Substituting this pressure 
into the expression for F adv and setting v = vk, we find that F adv ~ 13x~' 1 x~ l / 2 {g c c/ k)£ c ^ 
where g c is the gravitational acceleration at the top of the disk core. The total radiation 
flux passing through the halo, F tot , is equal to (g c c/ k)£ c . Under these assumptions we would 
therefore conclude that most of the energy — at least in the inner parts of the disk, say, 
near x±/ 2 ~ 23 where most of the energy is generated — is advected. If the advective flux 
dominates then the radiation becomes trapped in the outflow and the mass loss is maximized. 
The mass loss rate per unit disk area is then m w ~ F tot /v 2 K , which is just enough to prevent 
the total luminosity from exceeding the Eddington limit by a substantial factor. 

However, this argument appears not to be self-consistent. According to our prescrip- 
tion for photon bubble transport, the diffusive component of the flux is given by F diS = 
— (c/K,)(Vp)/p)£ ~ — 0.3xczVp. Unless x is small, at x 3> 1 the diffusive flux should domi- 
nate over the advective flux if z\Vp\ ~ 0{p). The reason we found that the advective flux 
could dominate in the halo solution is that, according to this solution, z\Vp\ <C 0{p). But 
as soon as the advective flux becomes appreciable, the acceleration increases and one can 
show from the equation of motion that this ordering no longer applies. We conclude that 
the flow cannot become strongly advection-dominated. 

Ignoring advection allows us to construct a simple, quasi-one-dimensional model for the 
flow dynamics. Defining y = z/r, we mimic the transition from vertical to radial flow by 
defining the conserved mass flux per unit disk area as m w = p(y)v(y)(l + y 2 ). By the same 
token, we write the radiation flux as -F to t ~ -^diff ~ (dcC/ '^)^c(l + y 2 ) 1 and the effective 
gravity as (v 2 K /r)y/(l +y 2 ) 3 ^ 2 . If we further define w = v/vk, Q = 8 c £ c VK/(0-3rh w Kr), and 
denote differentiation with respect to y by a prime, we obtain the equation of motion 

w v 

ww' = Q \ , (48) 

y (i + y 2 Y /2 

which is valid where £ > 1. This equation reproduces the structure of the static halo in the 
subsonic zone at y 1, with p oc z~ 2 implying v cx z 2 . 

Although this equation is tricky to solve exactly, because of the singularity, it is straight- 
forward to show that sensible solutions exist only for Q ~ 0(1). For Q ^> 1 (small rh w ), 
the inertial term on the left-hand side of eq. (48) is negligible for all y < 1. The flow never 
reaches escape speed and decelerates at large y. For Q 1 (large m w ), the inertial terms 
become appreciable at y ~ Q 1, but because the porosity due to photon bubbles (as 
measured by £) decreases with decreasing density, the flow never receives enough momentum 
to accelerate it to escape speed — again the flow decelerates. Only for Q ~ 0(1) does the 
flow receive the right momentum to accelerate it to ~ vk at y ~ 1, at which point it can 
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coast (with possibly some additional mild acceleration) to large r. 
The mass loss rate per unit disk area is then given by 

3 ^M E mD 

™- ~ < 49 > 

Mass loss begins at an outer radius x where the accretion rate m first exceeds the Eddington 
limit, m ~ 2xo/3. For the coefficients obtained above (taking into account mass loss 
from both sides of the disk and assuming that rh does not exceed the limit imposed by 
photon bubble transport), mass flux is depleted by the wind according to the equation 
drh/dx ~ l.5Q(T>/x 3 / 2 )rh. For x- m = 6 the solution is 



rh(x) 
m 



> exp [Q [3.7 {x- 1 - x, 1 ) - 3 (x- 1 ' 2 - x 1/2 )] } . (50) 



For Q — 1 and x < 1000, more than 50% of the accreting matter reaches the center, implying 
that highly super-Eddington luminosities are possible without destroying the accretion disk. 
This result should be regarded strictly as illustrative, because the coefficient in eq. (49) is 
very approximate and could be off by a factor of several in either direction. Moreover, the 
fact the Q ~ 1 suggests that F a ^ v is at least comparable to the diffusive flux. If F^v is a 
sizable fraction of F tot then the diffusive flux is reduced by a factor 1 — F a< x v /F tot . The same 
factor would multiply the first term on the right-hand side of eq. (48). Acceptable values 
of Q would be increased by the reciprocal of this factor, and the mass loss rate would be 
correspondingly diminished. This argument should be valid provided that the total radiative 
luminosity is approximately conserved through the outflow, i.e., if only a small fraction of 
the radiation flux is converted to kinetic energy. 

Even if the mass loss is not catastrophic, the fact that it is a substantial fraction of 
the accretion rate means that the total kinetic luminosity of the wind. L w , is likely to be 
only a few times smaller than the radiative luminosity of the disk. For super-Eddington 
luminosities this can represent a large flux of energy, which could dramatically affect the 
environment. 



3.4. Disk Temperature 

Normal accretion disks have color temperatures that exceed their effective temperatures 
because the radiation pressure at the "true" photosphere, where the escaping radiation is 
thermalized, is larger than the radiation pressure at the scattering photosphere. Typical 
color corrections, / co i = T color /T eff , lie in the range ~ 1.5 — 2 (Shimura & Takahara 1995; 



-26- 



Davis et al. 2005), implying that thermalization occurs at a scattering optical depth ~ 5 — 20. 
How does photon bubble transport affect disk color? 

A detailed treatment of this issue is beyond the scope of this paper, but we are able 
to make some comments about the expected trends. First, we note that the absorption 
opacity used to estimate thermalization effects should be evaluated at a density close to the 
maximum density found in the photon bubbles, rj + po, where po is the mean density. Although 
the densest gas occupies only a fraction ~ r/^ 1 of the volume, the absorption opacity /t a b s is 
proportional to r] + . Therefore the effective opacity is k cS ~ [77+ 

K cs K abs(po)] 1 ^ 2 where K es is 

the electron scattering opacity. The absorption opacity (presumably a mixture of free-free 
and bound-free opacities) generally decreases with radiation pressure in LTE — for free- 
free absorption we have /% oc p~ 7//8 , with a somewhat slower mean decline if bound-free 
transitions are important. 

First consider the disk core. Compared to a standard disk at radius x with the same 
m and a, the photon bubble-dominated disk has both higher pressure and higher density, 
p c oc £ c and p c oc t\. We also have i] + oc £ c , therefore the effective opacity is higher than in 
the equivalent standard disk by a factor i 25 ^ 16 . The effective optical depth for thermalization 
scales as T e s oc pK, c gH c oc 1 2 c k c r oc . Thus, thermalization should be much more effective 
in photon bubble-dominated disk cores than in standard accretion disks, mainly because the 
typical densities are much higher. 

Because the density and pressure drop by roughly the same factor between the core 
and halo, effective opacities (oc p 1 / 2 ^ 7 / 16 ) are unlikely to change by a large factor across 
the transition region. However, the large density drop (by a factor ~ 10a5 c / x 1 ^ 2 ) implies a 
substantial decrease in the effective optical depth. Whether the radiation is still thermalized 
in the lower parts of the halo therefore depends on details. Toward the upper halo, the 
pressure is roughly constant, p oc z~ 2 , and according to our estimates for case (c) photon 
bubble transport, 77+ oc z. Therefore, the effective optical depth decreases oc z~ 3 ^ 2 . Ther- 
malization will certainly fail well below the top of the halo, where the scattering optical 
depth approaches unity. These general arguments should also apply to the wind, in the 
super-Eddington case. If thermalization is still effective in the lower parts of the halo, the 
color correction should be close to unity, i.e., photon bubble-dominated disks would be cooler 
than their standard counterparts of the same luminosity, by a factor of ~ 1.5 — 2. 

The higher luminosities of super-Eddington disks, of course, would offset a smaller color 
correction. Roughly, we expect T color oc / co iM _1 / 4 (L/L E ) 1 / 4 . Therefore, a 1OM black hole 
radiating at lOL^ with no color correction would still appear hotter than a 100M Q black 
hole radiating at the Eddington limit with / co i = 2, but only by a factor of 1.6. If some 
thermalization occurred further out in the wind, or if radiation trapping were important in 



-27- 



the outflow, then the color temperature of the low mass hole could be smaller. 

4. Summary and Discussion 

4.1. Principal Results 

We have attempted to refine and flesh out our earlier treatment (B02) of the possible 
role of photon bubbles in accretion disks. The speculative nature of this work cannot be 
overemphasized, especially given the fact that it is unknown whether highly nonlinear pho- 
ton bubbles develop in the presence of hydromagnetic turbulence generated by MRI. We are 
motivated to take an optimistic view by the fact that photon bubbles develop in uniform 
fields on timescales much shorter than the timescales that characterize MRI (Gammie 1998; 
B02; Turner et al. 2005). Moreover, the fact that MRI appears to predict a predominantly 
horizontal magnetic field suggests that field coherence may be adequate to support photon 
bubbles of moderately long wavelength (and therefore of large density contrast). Never- 
theless, the time-dependence of the field, as well as its small-scale structure, could very 
well prove fatal to the premises upon which this paper is based. At least, it is hoped that 
the arguments presented in this paper will help to frame future radiation hydrodynamical 
simulations, which will settle these questions. 

4-1.1. Photon Bubble Theory 

Our analysis is based on an analytic model of photon bubbles as a periodic train of 
planar, isothermal shocks mediated by gas pressure (B01). The magnetic field is sufficiently 
stiff that it enforces one-dimensional motion. The general validity of this model has been 
verified by numerical simulations (Turner et al. 2005). However, photon bubbles in accretion 
disk cores occupy a part of parameter space that had not been analyzed previously. 

First, the optical depth across low density regions of the photon bubble, as measured 
by the parameter To — the optical depth for gas at the mean density, across the gas pres- 
sure scale height — is predicted to be small. The usual diffusion law relating flux to the 
radiation pressure gradient and the density has to be modified, since the flux cannot adjust 
to changes in density across distances shorter than a photon mean free path. We therefore 
generalized the closure to the radiation hydrodynamics equations described by Ruszkowski 
& Begelman (2003), to take account of fluid motions (see Appendix A). This generalization 
greatly complicated the equation describing the photon bubble structure (§ 2), changing it 
from a first- to a third-order nonlinear differential equation. Nevertheless, the generic flow 
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pattern, with a critical point at some characteristic flow speed, is preserved, and we were 
able to approximate the principal features of photon bubbles at finite and low optical depths. 
Most importantly, we found a tendency for the low-density parts of the photon bubble to 
"fill in," so that the optical depth across any scale length never becomes much smaller than 
one. Another new behavior, which appears as long a radiation damping is sufficiently weak, 
is an increase in the photon bubble phase speed to values well above the gas sound speed. 

Second, radiation damping is predicted to be important in accretion disks. A damping 
parameter £ — essentially the ratio of the radiation diffusion time to the sound crossing time 
in the gas — was included in the B01 models, which apply (with some minor corrections: see 
Turner et al. 2005) for To 1. We extended the analysis to include the effects of damping 
at finite and small optical depth. Although we were not able to provide a rigorous proof, 
it appears that photon bubbles with large density contrasts can form only when £ < 1 or, 
more precisely, when £ is less than some value which depends on the orientation of the shock 
fronts and the magnetic field but which is probably ~ 0(1) in practice. As in the case of 
photon bubbles at large r, the "Eddington enhancement factor" £ — which is the ratio of 
the flux to the local Eddington flux, gcj k, and not to the global Eddington limit — is given 
essentially by the ratio of the mean density to the minimum density, as first pointed out by 
Shaviv (1998). 



4-1.2. Accretion Disk Theory 

We applied our simple photon bubble theory to a one- zone a— model for the equato- 
rial region of the disk around a black hole. The treatment resembles standard models for 
radiation pressure-supported disks except for the inclusion of £ as a free parameter. Where 
photon bubbles do not exist, £ = 1. In the absence of photon bubbles, £ increases with r, im- 
plying that damping prevents the formation of photon bubbles beyond a certain radius. For 
accretion disks less luminous than about 0.3(a/0.03)~ 9 ^ 8 (M/M Q )~ 1 ^ 8 L E (for an accretion 
efficiency of 0.1), £ > 1 at all radii and photon bubbles cannot form. 

For more luminous disks, the onset of radiative transport by photon bubbles, at the 
radius where £ first drops below ~ 1, causes £ to increase with decreasing r. In a steady 
state with constant rh, unchecked energy loss via photon bubbles would increase £ to the 
point where photon bubble formation shuts off again. We therefore deduce that there is a 
feedback mechanism, which should regulate the amount of energy transported by photon 
bubbles so that £ ~ 1. The photon bubble density contrasts required to maintain such an 
equilibrium are small compared to the maximum values permitted by photon bubble theory, 
lending credence to the conjecture that such contrasts could be attained under turbulent 
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conditions. 

This feedback mechanism, which should operate on a dynamical timescale, makes spe- 
cific predictions about the thermal and viscous stability of disks dominated by photon bubble 
transport. Under the assumptions of an a— model viscosity, we find that the steady state so- 
lution with £ ~ 1 does not cure the well-known thermal and viscous instabilities that plague 
radiation pressure-dominated disks. An alternative outcome, in view of this instability, is 
that there is no steady accretion. Instead, the inner disk might cycle through a sequence of 
states in which the surface density builds up at small £, where the disk is thin and viscous 
transport inefficient, and is subsequently depleted by accretion or evaporation of a thickened 
disk, without photon bubbles, at large £. Recall from § 3.1.1 that both thermal and viscous 
instabilities can be quenched by photon bubble transport if £ 1 and the photon bubbles 
are able to approach the maximum permitted density contrast. 

Steady state disk cores with £ ~ 1 can liberate super-Eddington luminosities via photon 
bubble transport. This result is qualitatively and quantitatively similar to the conclusion we 
reached in B02, although the details of the photon bubble transport model, especially the 
feedback effect, are different. The predicted maximum luminosity depends weakly on mass, 
with L max ~ 2O(«/O.O3)(M/M ) 1 / 9 L £ . 

Super-Eddington luminosities necessarily drive winds. As the disk density declines, its 
ability to support photon bubbles with large density contrasts decreases. As £ decreases, 
the mean force exerted by the escaping radiation increases, driving the wind. In particular, 
if the advected flux Apv exceeds the diffusive flux, then the mass loss will probably deci- 
mate the accretion disk, decreasing m at such a rate (m oc r) that the luminosity cannot 
significantly exceed Le- Our approximate analysis suggests that the diffusive flux probably 
prevails, allowing most of the accreted matter to reach the inner disk. Although this argu- 
ment is tentative and will require detailed confirmation, it suggests that super-Eddington 
luminosities need not be compromised by mass loss. Mass loss could also be reduced if gas 
becomes trapped on closed flux loops anchored in the disk (B02; Socrates & Davis 2005). 

4.2. Caveats 

We have already emphasized the overwhelming caveat of this paper: that the incoherence 
and time-dependence of the magnetic field might prevent the formation of photon bubbles 
in accretion disks altogether. In this section we describe a number of other significant 
uncertainties related to the arguments presented. 

1. We used a highly approximate radiative transfer model to estimate the effects of fi- 
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nite optical depth on photon bubbles. The closure proposed in RB was tested against Monte 
Carlo calculations of radiative transfer through plane-parallel slabs. This test provides a rea- 
sonable facsimile of conditions in photon bubbles, and gives reasonable agreement (although 
the agreement is somewhat worse after the correction of an error in RB; see Appendix A) 
However, for the present problem we were required to generalize this closure to include the 
effects of motion. No comparable tests of been done of the generalized model. 

2. We did not solve the full, third-order nonlinear differential equation for photon bub- 
ble structure in the presence of finite optical depths. The full equation is two-point boundary 
value problem with a critical point and eigenvalue, in which the boundary conditions depend 
on each other and on the eigenvalue. To obtain a quick estimate of the functional dependence 
of the solutions on parameters, we devised an approximate method to find asymptotic solu- 
tions close to solutions that cross the critical point. We have no way of knowing how accurate 
this scheme is, although it gives reasonable (though not perfect) results in the high-r limit. 
We were also unable to prove that no high-contrast solutions exist for large £, although the 
approximate method suggests that this is the case. 

3. We adopted a simple, one-zone model for the disk core structure, assuming that 
the density scale height coincides with the thickness of the dissipation layer at each radius. 
However, simulations suggest that much of the dissipation may occur outside the density core 
(Miller & Stone 2000; Turner 2004; Hirose et al. 2005). If the bulk of the dissipation occurs 
at high optical depth, but outside the density core, then £ will be systematically lower than 
our estimates indicate. Photon bubble transport could then be important at larger radii 
and larger accretion rates. In this sense, our one-zone model is a "worst-case scenario" for 
assessing the potential importance of photon bubble transport. On the other hand, if most 
of the dissipation occurs at such low optical depths that photon bubbles cannot form (e.g., 
Socrates & Davis 2005), then the issue of photon bubble transport is moot. 

4. A related issue concerns the likely absence of photon bubbles deep within the disk 
core. £ is a monotonically decreasing function of height, formally blowing up at the equator 
where the vertical gravity vanishes. This suggests that in a disk core with £ c ~ 1, photon 
bubbles could exist only in the surface layers. What happens if a sizable fraction of the 
dissipation occurs near the equator, where photon bubble transport is absent? If this energy 
were carried to the photon bubble-dominated zone by radiative diffusion, it would puff up 
the disk core to the point that super-Eddington luminosities would be impossible. However, 
this is probably not what happens. The existence of the photon bubble-dominated zone acts 
like a drop in opacity near the surface, which causes the specific entropy to decrease with 
height in the disk interior and can drive strong convection down to the equatorial zone. We 
demonstrate this effect explicitly in Appendix B, by constructing and solving a toy model. 
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This model shows that the ratio of convective to diffusive flux is determined by effective 
porosity of the disk near the surface; for photon bubble transport the ratio is ~ t f. 
As long as the viscosity parameter a is C 1, the required convective motions are highly 
subsonic and the convection can be efficient (see, e.g., Blaes & Socrates 2003; Turner et 
al. 2005). There is no problem with convection transporting a super-Eddington flux. The 
damping parameter £ is also affected by the orientation of the shock fronts, represented by 
the parameter \i with < fj, < 1. We arbitrarily set /i to a fixed value associated with the 
fastest growing linear photon bubble modes in the high-r limit. Although this assumption 
seems consistent with large- r simulations at low £ (Turner et al. 2005), it is possible that 
the dominant orientation would be affected by damping, or that /i would adjust to local 
conditions so that £ remains small over a much wider range of conditions than we have 
assumed. 

5. Likewise, we adopted a very simple prescription for the effects of photon bubbles in the 
halo and wind, assuming that £ would be close to the maximum permitted value where £ < f , 
but would rapidly decrease with £ for £ > 1. This was the basis for the proposed feedback 
mechanism, but it ignores the dependence of damping effects on the relative orientation of the 
magnetic field and the photon bubble shock fronts. Also, we adopted a simple prescription 
for determining the maximum £, based on the maximum allowed photon bubble wavelength 
or the maximum density contrast permitted by the magnetic field strength. 

6. Finally, we modeled only the density inhomogeneities due to photon bubble shock 
trains propagating in a static magnetic field, under the influence of the vertical component 
of gravity. There are other ways of creating inhomogeneities in a radiation-dominated ac- 
cretion disk, and since enhanced radiation transport is a generic characteristic of highly 
inhomogeneous atmospheres (Shaviv 1998), we may have underestimated the disk porosity. 
In particular, the accelerations associated with MRI would create transient zones of artificial 
gravity that could lead to the formation of photon bubbles close to the disk equator. Ac- 
celerations associated with magnetic tension on curved field lines — also driven by MRI - 
would create density inhomogeneities via a mechanism akin to the Parker (1966) instability. 
Inhomogeneities with density contrasts of order the ratio of the magnetic pressure to the 
mean gas pressure are seen in radiation-MHD shearing-box simulations (Turner, Stone, & 
Sano 2002; Turner et al. 2003; Turner 2004); the extent to which these result from photon 
bubbles remains to be determined. 
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4.3. Astrophysical Consequences 

In B02, we discussed several applications of photon bubble-dominated accretion disks, 
most notably, ultraluminous X-ray sources (ULXs: Fabbiano 1989; Colbert & Mushotzky 
1999; Zezas & Fabbiano 2002) and rapidly accreting supermassive black holes. Here we 
briefly amplify on these suggestions, based on our new results. 

It remains controversial whether ULXs are intermediate-mass black holes radiating at 
sub-Eddington luminosities, or stellar-mass black holes radiating at a few times Le (or 
whether there are two classes: Watarai et al. 2005). Our new analysis suggests that several 
aspects of photon bubble-dominated disks may make it difficult to distinguish these cases. 
At a given luminosity, the higher color temperatures associated with smaller disks around 
lower mass black holes would tend to make super-Eddington disks appear hotter. Some 
ULXs appear hotter than typical black-hole accretion disks (Makishima et al. 2000; Mizuno 
et al. 2001) while others appear cooler (Miller, Fabian, & Miller 2004). However, we saw in 
§ 3.4 that the color corrections of photon bubble-dominated disks could be smaller than those 
of standard disks. This would partially compensate for the smaller surface area, although 
the super-Eddington disks probably would still be warmer. Radiation trapping effects - 
although we expect them to be modest — would work in the same direction, with photon 
energies adiabatically shifted to lower values. Another possibility is that the accretion rate 
at large r exceeds the maximum allowed by photon bubble transport. Such disks could also 
produce super-Eddington luminosities, but with m decreasing toward smaller r according 
to eq. (43). The energy output would be much less centrally concentrated than in a mass- 
conserving disk, and the spectrum should be considerably cooler. In this case, the mass loss 
would be driven directly from the thick disk, without the intermediary halo. Radiation trap- 
ping could be much more important, which could reduce the color temperature further. We 
note that a similar model was proposed by Watarai, Mizuno & Mineshige (2001). Finally, we 
point out that the winds from super-Eddington disks would have scattering photospheres 
at radii several times larger than their points of origin. This could smear out variability 
at the high frequencies associated with QPOs from stellar-mass black holes, although lower 
frequencies would not be affected (Strohmayer & Mushotzky 2003). 

Super-Eddington accretion could lead to the very rapid growth of supermassive black 
holes at high redshifts. Accretion rates could exceed the Eddington value by 1-2 orders 
of magnitude (depending on black hole mass), leading to e-folding times shorter than the 
Salpeter time by 1-2 orders of magnitude, i.e., < 4 million years for seed black holes of 
stellar mass, decreasing to a few hundred thousand years for M > 10 6 M Q . Although we 
have argued that mass loss need not decimate a super-Eddington accretion disk, the wind 
would probably carry an appreciable fraction of the luminosity in kinetic energy, and would 
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therefore have a substantial effect on the surrounding protogalaxy. The release of so much 
energy on such a short timescale (shorter than the dynamical time of all but the innermost 
regions of the galaxy) would be like setting off a bomb in the middle of the protogalaxy, 
rather than a quasi-steady wind. 

I am indebted to Ari Socrates and an anonymous referee for comments that helped to 
improve the paper. This research was supported in part by NASA Astrophysical Theory 
Program grant NAG5-12035 and by the National Science Foundation under Grants AST- 
0307502 and PHY-990794. Part of this work was carried out at the Kavli Institute for 
Theoretical Physics at the University of California, Santa Barbara; I thank the members of 
KITP for their hospitality. 



A. Radiative Transfer in a Moving, Inhomogeneous Atmosphere 

In this appendix we generalize the radiative transfer model developed by Ruszkowski & 
Begelman (2003; RB) to handle motion in an optically thick atmosphere containing large, 
localized density gradients. We also correct an error in RB. The approximation is intended 
to apply where the overall optical depth through the atmosphere is large, but where the 
optical depth across a density gradient length scale may be small. 

We approximate the intensity in the frame comoving with the gas as 

7(x,ft) = J (x,6) + — ft-F(x), (Al) 
An 

where F is the local flux vector in the comoving frame, Cl is the directional unit vector, and 
x is the position vector. We assume that all the odd moments of 1$ vanish, i.e., J Qilodfl = 
J QiQjQklodfl = 0. To 0(v/c), the equation of radiative transfer for a scattering atmosphere 
may be written 



i di i - i - i r 

+ + A—tt ■ V(fi • v) = -/ + — / Idn, (A2) 

acDt a ac An J 



where a = pn is the scattering coefficient and DI / Dt is the Lagrangian time- derivative of 
the intensity (Mihalas & Mihalas 1984). The third term on the left-hand side comes from the 
Doppler shift and contains a factor 4 instead of 3 because / is the integrated over frequency. 
Taking the zeroth moment of eq. (A2) we recover the usual energy equation, eq. (7). The 
components of the first moment are 

1 DFj c,_ 4 3 



ac Dt 



- -(V • T)< - — — / Qi(Q ■ F)Cl ■ V(Cl ■ v)dtt (A3) 
a acAir J 
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where T is the radiation stress tensor, the components of which are 



Tij = - I QifljlodQ. (A4) 



c 



The closure relation for in terms of Fj and J = j- J /c?0 can be obtained by calculating 
the second moment of eq. (A2), assuming the form of the intensity given by eq. (Al). This 
leads to the equation for the radiation stress tensor: 

Tij = %■ + -1^ - A f q.q ^ n . V (n • V )dn, (as) 

3 5crc cte, J ac Dt ac J 

where u = 4nJ/c is the energy density and the radiation pressure is u/3. The Eddington 
factor is 1/3 because the overall optical depth of the atmosphere is large. Using eq. (A5) in 
eq. (A3) (with the assumption that Jo is isotropic) then gives a set of second-order partial 
differential equations for the components of F, which we use in § 2. 

Note that the second term on the right-hand side of eq. (A5) differs from RB eq. (6), 
which we believe to be in error. This error propagated into RB eq. (7), and thence into 
the slab models used for comparison with Monte Carlo simulations. The function a 2 (fi) in 
these models should simply be the factor 5. We note that analytic fits using the corrected 
formulae are still reasonable, although fits using 10 instead of 5 appear to be much better. 
We use 5 in the expressions presented in this paper. 

The last term on the right-hand side of eq. (A5) represents radiation viscosity. If the 
optical depth across any flow structure, At, satisfies At <C c/v, then we can neglect this 
term as well as DF/Dt and the final term in eq. (A3), the DT / Dt term in eq. (A5), and the 
"radiation trapping" (Dp/Dt) term in eq. (7). This approximation is relevant in accretion 
disks, hence we use 

V • F + ApV • v = (A6) 
F = --(V-T) (A7) 
for the energy and radiative transfer equation, respectively, where 



B. Convective Transport in a Photon Bubble-Dominated Disk: A Toy Model 



Generic models of radiation pressure-supported disks tend to be convectively unstable. 
This tendency is exacerbated if there is a layer at high z where photon bubbles form, since 
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this has the same effect as introducing an opacity that decreases strongly with height. Since 
the radiation pressure gradient is determined by hydrostatic equilibrium, and the divergence 
of the radiation flux is fixed by the dissipation rate, a sudden decrease in the opacity leads 
to a strong entropy inversion. 

To estimate the magnitude of this effect, consider the vertical structure of a radia- 
tion pressure-dominated slab in the presence of gravity g{z) = —z and radiative diffusivity 
c/k — x{ z )- Since we will only be interested in ratios of fluxes, we will use convenient nondi- 
mensional parameters. The presence of photon bubbles at large z but not near the equator 
corresponds to dx/dz > 0. Assume that the emissivity is proportional to density and that 
the total flux can be decomposed into a diffusive part and a convective part, F = Fd + F c . 
The vertical structure equations are: 

p' = -pz (Bl) 
F' = F d + F^ = p (B2) 

F d = -X-. (B3) 
P 

where a prime denotes differentiation with respect to z. Equations (Bf ) and (B3) then imply 
Fd = X z - Any region with x = const, is convectively unstable, since F c = implies that p is 
constant while p must decrease monotonically with height. Moreover, any sharp increase of 
X with height would exacerbate the instability, since it would imply a sheet of high density 
suspended in the atmosphere (according to eq. [B2]), i.e., a density inversion. (Strictly 
speaking, the upper part of the dense layer might be locally stable, if p decreases steeply 
enough, but most of the atmosphere would be unstable.) Therefore, it seems reasonable to 
suppose that the model atmosphere is convective at all z. If convection is efficient, then we 
can adopt the equation of state p = Ap 4 ^ 3 and the structure of the atmosphere is independent 
of x(z): 

p = p (l - x 2 f; p = p (l-x 2 )\ (B4) 

where x = z/z and z = (8Ap^ 3 y/ 2 is the height of the atmosphere. The total column 
density of the atmosphere is 

S = Po z I (1 - x 2 fdx = 0A57p z . (B5) 
Jo 

Suppose that x( z ) satisfies the boundary conditions x(0) = 1 and x( z o) — Xo ^ 1- 
Integrating eq. (B2) we have F c (z ) = S — Xo z o- But we also expect that all of the flux 
near the upper boundary should be carried by radiative diffusion, not convection, so that 
F c {zo) = 0. This implies p = 2.19x - 
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Now, consider the fluxes carried by the atmosphere near the equator, i.e., at x <C 1. 
Since x ~ 1 we have F d m z x while F c (p — 1)^0^ — (2-19xo — l)z x. Therefore the ratio 
of convective to diffusive flux near the equator is 

^(x< l)w 2.19x0-1- (B6) 

If xo represents the porosity associated with photon bubbles near the top of the disk, then 
Xo ~ ( 3> 1, where I is the local Eddington enhancement factor. The convective flux is 
therefore much larger than the diffusive flux, and is super-Eddington with respect to the 
local gravity at z <C z . This simple model shows that the development of photon bubbles 
near the photosphere can drive very strong convection near the equator. Interestingly, even 
when photon bubbles are unimportant, Xo — 1> the convective flux near the equator is 
slightly larger than the diffusive flux in this toy model, F c /F d (x <C 1) ~ 1.2. 
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